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This course has a dual purpose. First we review the successes of the weak-coupling BCS theory 
in describing new classes of superconductors discovered since 1979. They include the heavy-fermion 
superconductors, organic superconductors, high- Tc cuprate superconductors, Sr2Ru04 etc. Second, 
we present the quasiclassical approximation introduced by Volovik, which we extend to describe the 
thermodynamics and the thermal conductivity of the vortex state in nodal superconductors. This 
approach provides the most powerful tool in identifying the symmetry of the energy gap function 
A(k) in these new superconductors. 

I. INTRODUCTION 

Question: What is the difference between a Fermi liquid and a non-Fermi liquid? 
Answer: The difference is the same as the one between bananas and non-bananas. 

Boris Altschuler (2001) 

Unconventional or "nodal" superconductors appeared on the scene in 1979, when the heavy-fermion superconductor 
CeCu2Si2 and the organic superconducting Bechgaard salts (TMTSF)2PFg were discovered. Since then, many more 
heavy-fermion superconductors Q and organic superconductors have been synthesized. This development was 
followed by the epoch making discovery of high- Tc cuprate superconductor La2_a;Ba2;Cu04 by Bednorz and Miiller 
I^J with the superconducting transition temperature — 35K in 1986. Within a few years new classes of high-T^ 
cuprates emerged, including YBa2Cu306+5, La2-a;Sr2:Cu04, Bi2Sr2Cai_a;Ya;Cu208+5 (Bi2212); and HgBa2CaCu206 
with Tc = 145K. The subsequent enthusiasm and confusion are well documented in an early review by Enz j^. 
Confusion? Yes, initially it was thought that Landau's Fermi liquid theory Q and the BCS theory Q were no longer 
applicable. 

Among many proposals one of the most influential were Anderson's dogmas "T], which can be summarized as follows: 

a. The action takes place in the Cu02 plane common to all high-Tc cuprate superconductors. 

b. The undoped state is a Mott insulator with antiferromagnetic (AF) order. Upon doping superconductivity 
appears. Therefore the simplest Hamiltonian is the two-dimensional (2D) one-band Hubbard model: 

H = -tY,i4c.CJc. + h.c) + uY,n^'trl,^, (1) 

{i,j)a i 

where {i, j) connects the nearest neighbors in the 2D square lattice. The and Cia are creation and annihilation 
operator for the holes at the site i with spin a and ~ c^Q^. 

c. As a possible ground state of Eq. Anderson proposed the resonating valence bond (RVB) state: 

*=n(i-'^^)iBcs), (2) 

i 

where |BCS) is the BCS state for s-wave superconductors and Y\i{^ ~ '^i) with di = ni-^nn is called the 
Gutzwiller operator. Y\^{1 — di) annihilates all doubly occupied states. 

In spite of tremendous efforts spent on both Eqs . and 10) it has been difficult to find solutions in two dimensions. 
On the other hand the ID version of Eq. Q is now completely understood In the m eantime the perturbative 

analyses based on Eq. ^ predict BCS d-wave superconductivity in high- Tc cuprates |lCillll2il j |. 



FIG. 1: Local density of states around a single vortex for d-wave (left) and s-wave superconductivity. 



High-quality single crystals of YBCO, LSCO and thin films of Bi2212 became available around 1992. The d-wave 
superconductivity in these high- Tc cuprates was established in 1994. A mong many experiments, angle-resolved pho- 
toemission spectroscopy (ARPES) [Tjl and Josephson interferometry p!5llla | played a crucial role. Around this time 
several theoretical groups started to analyze the physical properties of d-wave superconductivity within the BCS 
framework jl7i.l8.19..20j . In 1993 Patrick Lee demonstrated the universal heat conduction in d-wave supercon- 
ductors. Furthermore the thermal conductivity was shown to increase with increasing impurity scattering [13. This 
counterintuitive behavior was observed in the Zn-substituted YBCO in ^2^. The electronic contribution to the ther- 
mal conductivity in d-wave superconductors is proportional to T at low temperatures (i.e. T ^ A where A is the 
maximum value of the energy gap). Here we assume A(k) — Acos(2(/)) and (j> is the angle k in the a — 6 plane makes 
from the a axis. Then in the limit of no impurity scattering (i.e. F — > where F is the quasiparticle scattering rate 
in the normal state) the thermal conductivity is given by 



^ - (3) 

where kqo/T = limr^o f^/T, and v/v2 = Ep/^ and n is the quasiparticle density or the hole density. The velocities 
V and V2 are defined from the quasiparticle energy at the Dirac cone 



i?k - ^vH^~kFr+v^,ki, (4) 

where v is the Fermi velocity, fc|| is the radial component of the wave vector and k± the component perpendicular to 
the Fermi surface. 

Later we shall derive Eq.(3) in Section IV. The thermal conductivity in single crystals of optimally doped YBCO 
and Bi2212 below T = IK was measured by May Chiao et al ^23,24]. Making use of Eq. Q they found A/Ep = 1/10 
and 1/14 for Bi2212 and YBCO respectively. These remarkable ratios imply: 

a. High-Tc cuprates are described by the BCS theory of d-wave superconductivity. They are far away from the 
Bose-Einstein condensation limit which requires A ~ Ep [25| . 

b. According to the Ginzburg criterion, fluctuation effects are of the order A/Ep ~ {ppO^^- other wor ds they 
should be at most 10%. This appears to exclude large phase fluctuations and stripe phases discussed in 1291231 . 

c. For A/Ep — 1/1 there are hundreds of quasiparticle bound states around the core of a single vortex in d-wave 
superconductors |28l29| . The radial (r) dependence of the local quasiparticle density of states is very similar 
to the one obtained for s-wave superconductivity [s^. We show in Fig. 1 the local density of states around a 
single vortex for a d-wave and s-wave superconductor. These are well-known bound states first discovered by 
Caroli, de-Gennes and Matricon [3lll33 | for s-wave superconductors. 

In earlier works '33 34 35 36] it was asserted that there would be no bound state around a single vortex in d-wave 
superconductors. However in these works it was assumed A ~ Ep in order to facilitate the numerical analysis 
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FIG. 2: The phase diagram of hole-doped high-Tc cuprate superconductors. 

based on the lattice version of the Bogohubov de Gennes equation. This assumption (A ~ Ep) eUminates the 
bound states in these numerical analyses. If we limit ourselves to the experiments on single crystals of high-T^ 
cuprate superconductors, we find that most of the behaviors observed are consistent with the BCS theory of d- 
wave superconductivi tv l37l . Also it i s bette r to use the continuum version of the Bogoliubov de Gennes equation, 
which is proposed in |38j| and used in [28l29j |. More recently a similar analysis is extended for a vortex in an f-wave 
superconductor [3^ . 

From a theoretical point of view the universality of the Landau Fermi liquid in 2D systems was demonstrated within 
the renormalization group analysis (4Qii4L»4^] . The quasiparticles in the normal state of high-Tc cuprates appear to 
be a Fermi liquid state. Furthermore the quasiparticles in d-wave superconductors are in a BCS-Fermi liquid state 
with the quasiparticle energy 

= ^v^{k\\ - kp)^ + A^cos^{2<j>). (5) 

In the vicinity of the Dirac cone Eq. ||SJ) reduces to Eq. Q). 

Another consequence of the renormalization group analysis is that the instability of the normal Fermi liquid is 
related to the infrared divergence in the particle-particle (and the hole-hole) channel or the particle-hole channel. The 
former results in conventional or unconventional superconductivity and the latter in conventional or unconventional 
density wave states. Therefore it is of great interest to look at the phase diagram of the high-Tc cuprates from this 
point of view. We show a schematic phase diagram in Fig. |2J 

Let us consider the hole-doped region. In the vicinity of a; = there is an antiferromagnetic (AF) insulating phase, 
which is a Mott insulator (MI). As the hole-doping x is increased, the AF order is rapidly suppressed around x = 3%. 
Then a superconducting region develops for 5% < a; < 25%. This is sometimes called "the superconducting dome" . 

Also in the underdoped region there is the pseudogap (PG) regime. The nature of the pseudogap is still hotly 
debated. There is eviden ce that it is a d-wave density wave (dDW) ,43.44.45.46i | with an energy gap A(k) = A cos(20). 
Earlie r propo sals |43l45j | have only considered the commensurate case with Z2 symmetry something similar to a flux 
phase |47ll4^ . 

On the contrary, we consider the incommensurate c?-wave density wave |44l46j| case where the condensate has U(l) 
symmetry as in conventional charge-density wave systems (CDW). Here the phase vortex is the most common 
topolo gical defect. Moreover the phase diagram suggests a coexistence region of dDW and d-wave superconductivity 
|43l45l |. Then the phenomenological gap introduced by Tallon and Loram 50] should be the energy gap associated 
with dDW. 

Recently Laughlin |5l| has pointed out that the wave function (2) is impractical, since the Gutzwiller operator has 
no inverse. Instead, he proposed to analyze Eq. (|2Jl with a modified Gutzwiller operator which has an inverse. For 
example 

l[{'^-d^) ^ l[{l-ad,), (6) 



FIG. 3: From top left: 2D /-wave in Sr2Ru04, d^2_y2-wave in CeCoIns and in fi:-(ET)2Cu(NCS)2, (s + p)-wave in YNi2B2C, 
{p + /i)-wave in the A phase of PrOs4Sbi2, (p + /i)-wave in the B phase of PrOs4Sbi2. 

with a < 1. He called this "gossamer superconductivity". "Gossamer" means filmy cobweb or something light, fragile 
but strong. If we look at this new wave function with Eq. © , we realize that this is an example of competing order 
parameters |52l53j . We shall come back to this question at the end of our course. 

In a seminal paper Volovik |54j has shown that the quasiparticle density of states in the vortex state of the d-wave 
superconductors is calculable within a semiclassical approximation. The predicted ^/H dependence of the specific 
heat has been confirmed experimentally in sing le crystals of YBCO !55"5^, LSCO '57^ and Sr2Ru04 [slsgj . This 
semiclassical approach has been extended in a variety of directions .60 61 62 63 64,65.] . When high-quality single 
crystals in the extremely clean limit (i.e. / ^ £, where / is the quasiparticle mean free path and ^ the superconducting 
coherence length) are available, the angle- dependent thermal conductivity in the vortex state provides unique access 
to the gap symmetry A(k). In the last few years Izawa et al. have determined A(k) in Sr2Ru04 66], CeCoIns |67l ]. 
K-(ET)2Cu(NCS)2 iHl, YNiaBzC ^ and PrOs4Sbi2 |7nl7l| . These gap functions are shown in Fig. El 

In the following, we first focus on the quasiparticle spectrum in a variety of nodal superconductors. Then the effect of 
impurity scattering and the universal heat conduction is briefly summarized. The quasiclassical approximation in the 
vortex state in nodal superconductors is the central part of this course. Also the properties of nodal superconductivity 
in YNi2B2C and PrOs4Sbi2 are briefly summarized. In the last chapter we discuss unconventional density wave and 
gossamer superconductivity, which indicate new directions to follow in new materials. 
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II. QUASIPARTICLE SPECTRUM IN NODAL SUPERCONDUCTORS 

Following the BCS paper we consider the effective Hamiltonian 

^ = X!^(^)"=kaCka + ^ H ^(k^kOcJ-a-C^k' -a'C-k-aCka, (7) 

with 

z;(k,k') = -(|/P)-V/(k)/(k'), (8) 

and 

= i- 1 dn\f{k)\' 3D (9) 

dxdc^\f{k)\' 2D (10) 
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depending on whether the system is 3D or quasi-2D. In the following, we consider a group of quasi-2D superconductors, 
whose quasiparticle density of states, thermodynamics etc. are identical 64]. We define quasi 2D systems by a 
cylindrical Fermi surface, as in high-T^ cuprates, Sr2Ru04, CeCoIus, k-(ET)2Cu(NCS)2 etc. Also we consider a 
group of nodal superconductors with /(k) = cos(2</)), sin(20) (d-wave), e^**^ cos(x) (/-wave), cosx (d-wave), sinx 
(p-wave), e***^ sin(x) (c?-wave), cos(2x) (gi-wave) etc. Here x = ckz- Then it can be readily shown (Exercise 1) that 
these superconductors have an identical quasiparticle density of states (DOS). 

Within the mean-field approximation, i.e. the BCS approximation, Eq. Q is transformed as 

^ - E KM^)P3 + A(k)pi)vi/k,. - E (11) 

k,a k 

The corresponding Nambu-Gor'kov Green function is given by 



G-i(k,L^)=^-e(k)p3-A(k)pi, (12) 

where the pi's are Pauli matrices operating on the Nambu spinor space. For simplicity we consider here only spin 
singlet pairing and / as a real function. 

Then the poles of the Green function Eq. (|12|l give the quasiparticle energy 

Lo = ±v/e'(k) + A2(k) ~ ±^v^m - kpf + vlkl. (13) 

The last expression is an approximation near the Dirac cone. 

From Eq. (|13|) the quasiparticle density of states is obtained as ^3 

/ \E\ \ { f l^l^(^) for \x\ < 1, 
giE)^Re{^=U===)^\ (14) 
\ v/S2_ A2|/|2/ 1^ 2.K{l/x) for > 1, 

where x = \E\/A and K{k) is the complete elliptic integral of the first kind. We show in Fig. 01 the quasiparticle 
density of states versus x, compared with the one for s-wave superconductor with a full energy gap. 
For small energies, the density of states can be expanded as 

g{E) ^ \E\/A (15) 



for \E\/A < 1. 
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FIG. 4: The quasiparticle density of states for nodal superconductors (solid line) and for s-wave superconductors (dashed line). 



The mean-field approximation also gives the gap equation 



= {\f\' 



/ \E\ \ E 

dERe ( , ' ' = ) tanh(— ), 



(16) 



which we have written in 2 alternative forms. Here A = vNq is the dimensionless coupling constant, LUn is the 
Matsubara frequency and Eq is the cut-off energy. Also the cOn sum in the first equation has to be cut off at ujn — Eq. 
Then within the weak coupling limit we find |74| 



A(0)/r, - 2.14, 



(17) 



A(i)/A(0) ^ Vl^, 

where t — T/Tc- In Fig. Owe show A(<)/A(0) versus t together with the approximate expression Eq.(18). 
As to the thermodynamics, it is convenient to start with the entropy S given by 



/.oo 

5 = -4iVo / dEg{E)[f\n,f+[l-f)\n{l-f)] 
Jo 

= ANo r dE g{E)[\n{l + e"'^^) + (3E{1 + e'^^)]. 
Jo 



(18) 



(19) 



Here iVo is the quasiparticle density of states at the Fermi surface in the normal state and / = (1 + e^^)^^ is the 
Fermi function. 

From S, the specific heat and the thermodynamic critical field are obtained by 



Cs/T 



dS_ 

5T' 



(20) 



and 



-^ij2(r) = K(r) - Fs{T) ^ f 'dT' (5„(r') - 5,(r')). 



(21) 




FIG. 5: The temperature dependence of the order parameter (dashed - approximate, solid - exact 



Using the density of states g(E), we obtain 



C, ^ 27C(3) T 
jsT 27r2 A 



for T < A, 



(22) 



and 



iJc(O) = V2vrA^oA(0). 
Similarly the superfluid density within the a — & plane is given by 



(23) 



Ps 



{T)/ps{0)^27:TJ2 



2 j-2 



In the limit T ^ this reduces to 



(24) 



Ps{T)/p,{0) = l-21n2 



A(0)- 



(25) 



In other words one can obtain A(0) from the T linear slope of the superfluid density. 

There is still a controversy as to the correct expression of the c axis superfluid density. The simplest assumptions 
[TSf give 



PsAT)/PsAO) = ^||^(/tanh(A//(2T))) 

9 (26) 



6 [a{0)J 

Finally, the spin susceptibility and the nuclear spin lattice relaxation rate are given by 



Xs/Xn^l-Ps{T)/ps{0), 



(27) 
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and 



" ' (28) 



E 
^2T' 
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We stress again that these expressions are not only valid for d-wave superconductors in the weak-coupling limit 
as in high-Tc cuprates, but for all nodal superconductors with /(k) given above. Therefore in order to explore the 
individual gap symmetry, we have to look at other properties. 

It is somewhat surprising that all the energy gaps of the identified nodal superconductors in the quasi-2D systems 
belong to the above class of /'s. 



Exercises 1. 

1.1 Within the weak-coupling theory calculate the jump in the specific heat at T — Tc 

a. for s-wave superconductivity. (Answer: AC/C„ — 12/(2C(3)) = 1.43) 

b. for nodal superconductivity with A(k) = A/(k). (Answer: AC/C„ = 12/(7C(3))(|/|'')/(|/|2)2) 

1.2 Evaluate the ratio A{0) /Tc within the weak-coupling theory 

a. for s-wave superconductivity. (Answer: A(0)/Tc = n/j = 1.76) 

b. for nodal superconductivity. (Answer: A{0)/Tc = (7r/7) exp(— (/2)~-'^(/2 In/))) 

1.3 Show that the quasiparticle density of states for a group of /'s discussed is the same as given by Eq. (|14|l . 



1.4 Express Ke{f'^/y'x'^ — f^) in terms of complete elliptic integrals. Answer: 

/ f2 \ f UK{^)^E{x)) forx<l, 
Re ( ■' ) = <^ (29) 
W^^-fy [ ^x{K{l/x)- E{l/x)) fora;>l. 

III. EFFECT OF IMPURITY SCATTERING 

In metals the presence of impurities or foreign atoms is unavoidable. Also, they provide the simplest agents of 
quasiparticle relaxation. Therefore the study of impurity scattering is crucial to understand quasiparticle transport 
such as electric conductivity and thermal conductivity. 

In the early sixties the effect of impurity scattering in s-wave superconductivity was systematically studied in 
|7fil77l7a^ . As is well known, nonmagnetic impurity scattering has little effect in s-wave superconductors. The 
superconducting transition temperature and the thermodynamics are almost unaffected. The most dominant effect 
is a reduction of the quasiparticle mean-free path, as in the normal state, and its consequence on the magnetic 
penetration depth. On the other hand, magnetic impurities have a profound effect on s-wave superconductivity. The 
superc onducting transition temperature is sharply reduced. Also in some cases gapless superconductivity is induced 
|7Sl79j . ' ^^.^ 

In contrast, the nonmagnetic impurities have profound effects on nodal superconductors as pointed out in |8Ci8ll83 |. 
In particular, resonant impurity scattering appears to be prevalent. The extreme limit is the unitary limit where a 
resonance occurs at E = 0. Such a model has been discussed for c?-wave superconductivity in high-Tc c uprates in 
^3 84] . The first self-consistent studies of impurity scattering in d-wave superconductors were performed in |l8ll9l20j . 

For simplicity we shall limit ourselves to the unitary limit. Also we assume that the impurity is point like, i.e. it has 
only s-wave scattering amplitude. Then the effect of impurity scattering can be incorporated by w — > w in Eq. H12|) 
where oj is the renormalized frequency given by 
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FIG. 6: The quasiparticle density of states in the presence of impurities, for several impurity concentrations. 



and X — uj I IS. and T — ^^(TrA^o) ^ is the quasi-particle scattering rate in the normal state. Then the quasiparticle 
density of states in the presence of impurities is given by 

The quasiparticle density of states in the presence of impurities is shown in Fig. |3 We note a rapid appearance of 
disorder-induced spectral weight at £^ = 0. Indeed g{0, T) is given by 

5(0, r) = -j£^K (^^L=] , (32) 



where Co is obtained from 



^^--x^^-M^= ■ (33) 



In the hmit T/A 0, both Eq. ^ and Eq. ^ reduce to 



.(o,r)^Coin(±). gm- 




(34) 



and 



C2ln(4/Co)^^. (35) 
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FIG. 7: A(0,r)/Aoo (dashed line), Tc/Tco (solid line) and g{Q,V)/NQ (dashed-dotted line) are shown as a function of F/Fc in 
the unitary limit. 



The gap equation in the presence of impurity scattering is given by 



where Xn = Wn/A and a3„ is the renormalized Matsubara frequency. Then in the Umit A — s- 0, we obtain the 
Abrikosov-Gor'kov formula: 



where Tj, (T^o) is the superconducting transition temperature in the presence (absence) of impurities. Here ^'(z) is 
the digamma function. 

Note that Eq. (|37l) is the same as for s-wave superconductors in the presence of magnetic impurities 78] . In nodal 
superconductors T is due to nonmagnetic impurities, whereas in s-wave superconductors F is associated with magnetic 
scattering which involves spin flipping. 

At T = OK Eq. ^ reduces to 



- In [^^) - 2 ^/^ ln(Co + ^| + P)^ ' ^ Jc '^'''''^^ ^^^^^^^ ^ ""'^^ ~ ^^^^ 

where E — i?(l/vT+~z?) and K = K{l/\/l + x^) and Co has already been defined in Eq. Ip!^. We show in Fig. [3 
Tc/Tco, A(0,r)/A(0,0,) and g{0) versus F/F^ where F^ = nTco/{2j) ~ 0.8819Teo. 

Surprisingly, A(0, r)/A(0, 0) follows very closely Tc/Tco- Also 5(0, F) increases very rapidly with F. This rapid 
increase in the DOS has been measured by the low temperature specific heat in doped LSCO [83. Note Cs/{"/sT) = 
g{0). As seen from FIG. IHlthe experimental data agree very well with the theoretical prediction. Especially the 
agreement is almost perfect in the vicinity of the optimally doped LSCO. This suggests strongly that the weak-coupling 
BCS theory for c?-wave superconductivity is adequate for LSCO. We have already mentioned that A(0)/Tc = 2.14 in 
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the weak-coupling limit. We can deduce this ratio for optimally doped LSCO, YBCO and Bi2212. These are 2.15, 
2.64 and 5.3 respectively. Indeed LSCO appears to be in the weak-coupling limit whereas YBCO may be in the 
intermediate regime. Since the large ratio 5.3 for Bi2212 is rather strange, we wonder if another order parameter is 




0^ I 

FIG. 8; Experimental values for T^{V) for doped LSCO 

hidden to make this large energy gap possible (see Section VIII). 

The planar superfluid density in the presence of impurities is given by 



For T — this reduces to 



PsiT,T) F T f^dx f E 

p. (0,0) ACo ^ Aico V K 



(40) 



This sharp decrease in superfluid density in Zn-substituted YBCO has been observed by the /i-SR experiments 8.6|- 
In Fig.Elwe show p,(T, F)/p,(0, 0) versus T/Tcq. 



Exercises 2. 

2.1. Derive Eq. 

2.2. Derive Eq. ^ and Eq. 



IV. UNIVERSAL HEAT CONDUCTION 



In 1993, Patrick Lee |22j pointed out that the thermal conductivity in d-wave superconductors is linear in T for 
T ^ A(0), and that it takes the universal value kqq/T — kgvn/{3Hv2) in the limit F ^ 0. A further study indicates 
that k/T increas es very rapidly with F, the quasiparticle scattering due to impurities This result was used by 
May Chiao et al |23l24j to extract A(0)/£'f of optimally doped Bi2212 and YBCO. By transforming the expression 
of the thermal conductivity given by Ambegaokar and Griffin 87], the thermal conductivity for T <^ A(0) is given by 
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A \ (C2 + /2)3/2 



!)3/2 / 

(41) 



2T 1 1 



where Cq and A — A(0,r) have been given in Eq. H33() and Eq. H38|l respectively and E{k) is the complete elliptic 
integral of the second kind. For later purposes, it is more convenient to normalize Kxx by k^, the normal state thermal 
conductivity, by F = Ec = 0.882Tco the critical scattering where superconductivity disappears. 

< = P^- (42) 
ol cm 



Then we can rewrite Eq. (|41|l as 



2F 1 / 1 \ 

/< = -^-==E - /i(F/F,). (43) 



/i(F/Fc) versus F/F^ is shown in Fig. 11111 

We note that Kxx increases monotonically with F. In other words, the thermal conductivity increases with the 
impurity scattering. This counterintuitive behavior is understood, if one realizes that the impurity scattering produces 
quasiparticles due to the pair-breaking effect . Indeed, the predicted F dependence of Kxx was verified in Zn-doped 
YBCO ll^l and more recently in Sr2Ru04 [sll- If we substitute in Eq. (|jT|l / = sin(2(/)), cosx, e^^'^cosx, sinx, 
e^^'^sinx or cos(2x), we will obtain the same result. The planar thermal conductivity can thus not discriminate 
between different nodal superconductors 'S!^. The result for the out-of-plane thermal conductivity is of more interest. 
We obtain 

'^../<-^/(l-cos(2x)) ^° 



'^"{Cl + m^l^l (44) 
= A(F/F,) 




T/TcO 



FIG. 9: The superfluid density versus temperature for various impurity concentrations. 
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FIG. 10: The functions /i(r/rc), hiT/Fc), Hr/Tc) 



for / — cos(20), sin(20) and cos(2x). But 



^ 4r, 1 

= ^2(r/r,) 



(^2 



(1 + C'o)E 



CI 



- K 



(45) 



for / = cosx, e^*"^ cosx and 



'^^^^'^^ " (ttTw) " (ttTw) ) (46) 

for / = sinx and e^'-^sinx- We show hiT/Tc), hi^/Tc) and h{^/'^c) versus r/r^ in FIG. CHI It is clear that 
Kzz for / = sm(j) and e^^'^sincf) in the limit F ^ 0. There is no universal heat conduction. Very recently 
thermal conductivity data for n^x and Kzz a.t T — 0.4K and H || z in UPd2Al3 was reported [89ll90l |. This is shown 
in Fig. ^3 Although the field dependence of kh is not the same as its F dependence, we conclude A(k) in UPd2Al3 
is more consistent with A(k) cos(2x) rather than A(k) ~ cosx- The latter has been proposed in |£1,] based on 
antiparamagnon exchange with Q = (0, 0, tt/c) [9^ . 

Exercises 3. 

3.1. Derive Eq. H41() for a variety of /'s. 

3.2. Calculate k^z for these /'s. 
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FIG. 11: Thermal conductivity data for UPd2Al3 



V. QUASICLASSICAL APPROXIMATION 



A. Vortex state 



In the presence of a magnetic field, all known nodal superconductors enter into the vortex state when the magnetic 
field exceeds Hci{T), the lower critical field. In this vortex state the quantized vortex lines form a regular two 
dimensional lattice in the equilibrium configuration. Also each vortex line carries a unit flux = hc/(2e) — 2.06 • 
10~^Gcm^. The marvelous theory of type II superconductivity was created by Abrikosov [93] in 1957, based on 
the phenomenological Ginzburg-Landau theory 94]. The microscopic foundation of GL theory was provided by 
Gor'kov [93 after the appearance of the BCS theory In s-wave superconductors the vortex lattice is usually 
hexagonal. But in d-wave superconductors the square vortex lattice is more stable for T/T^ < 0.9 |9ML]- 

Indeed 

the square vortex lattice was recently observed in LSCO by neutron scattering experiments~j98l] . Also unlike in 
s-wave superconductors, the quasiparticles dominate the low-temperature transport properties in the whole vortex 
state in nodal superconductors. The quasiclassical approximation then provides the most practical way to handle the 
quasiparticle transport in the vortex state for T ^ A. 

B. Quasiparticle spectrum in the vortex state 

In 1993 Volovik 54] showed how to calculate the quasiparticle DOS in the vortex state of d-wave superconduc- 
torSjWhen T <C A(0). Here we shall follow this procedure. First we note that the quasiparticle energy is shifted 
to i?k — V • q = i?k — k • in the presence of a superflow Here v, 2q and Vg are the quasiparticle velocity, the 
pair momentum and the superfluid velocity respectively. Also for a class of nodal superconductors we considered in 
Section 11, we obtain g{E) ^ \E\/A for \E\ < A. 

In the presence of a superflow this is generalized as 

g{E,ll} = A-'{\E-vq\), (47) 

or 

<?(0,H)^5(H) = A-i(lvql), (48) 

where (...) denotes the averages over the Fermi surface and over the vortex lattice. Here v • q is called the Doppler 
shift. 
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FIG. 12: Quasi-two-dimensional cylindrical Fermi surface 



Let us consider d-wave superconductors in a magnetic field H j| c. Then Eq. (|48|) becomes 



(49) 



Here d = 1/^/eH and we assumed for simplicity a square vortex lattice with lattice constant d. Also, in the present 
calculation we took q = 0/(2r), where r is the distance from the center of a vortex and a is the angle between v and 
q. We treated the average over the vortex lattice a la Wigner-Seitz. Finally in the earlier treatments 60 61 63 6^ the 
factor 7r~^ coming from {2tt)-^ J^"" d(t)S{cos{2(P)) was missing. From Eq. (|49|l the specific heat and other observables 
llOd are obtained as 



Cs/'fnT ^ g(H.), Ks/K„=g(H), 



Ps{0) 



gin). (50) 



We have already mentioned that the V H depen dence of the specific heat in the nodal superconductors was observed 
in YBCO [sslls^ . LSCO [E^l and Sr4Ru04 |58ll59| |. However, in the earlier analysis [sslls^ a factor ~ 0.3 was missing. 
We believe the tt"^ we find here accounts for this missing factor. In other words, the semiclassical result is not only 
qualitatively but also quantitatively accurate. It is very easy to work out (|v • q|) for other classes of A(k) as discussed 
in Section II. Then it is easily seen that the configuration H || c cannot discriminate d-wave superconductors from 
other nodal superconductors. 



C. Extension to the arbitrary field orientation 



In order to get a handle on the gap symmetry of A(k), it is necessary to consider the case of arbitrary field 
orientations. For simplicity we limit ourselves to the case where H is in the a — b plane with an angle (f) from the 
a axis. Precisely this configuration is considered in |^. Unfortunately, however, a ra ther unrealistic Fermi surface 
was considered. Instead we consider the quasi-2D Fermi surface shown in Fig. 1121 [lOlj . In the present configuration 
the vortex loses circular symmetry around the vortex axis. Therefore we find it is very useful to introduce a scale 
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transformation as in |102| to make a circular vortex. For / = cos(20) we find |6; 

2 



(7(H) = — wv eH ( y sin x + sin' 



2 I •„2^ 

(51) 



2 , 

~ -wV^(0.955 + 0.0285cos(4(/))), 

TT 

where v = ^/v^ and Wc is the Fermi velocity parallel to the c axis and (j) is the angle H makes from the a axis. We 
thus find a 3% fourfold term (~ cos(4(/))) with a <.3% angular variation of g(H) The (p dependence of the specific 
heat in YBCO was studied by Wang et al |103| |. They could not find this small fourfold term which is within their 
experimental error. 

When A(k) has a horizontal node as in Sr2Ru04 5(H) does not exhibit the (f> dependence. In particular for 
/ = e^'"^ cos X we obtain 



(|vq|> = ^vV^V2E{l/V2) = 1.216\iVeH (52) 

Therefore the specific heat of Sr2Ru04 by Deguchi et al |lfl4| appears somewhat puzzling. Their data for T > 0.5K is 
consistent with 2D /-wave superconductor. However, something different appears to happen be low T = 0.3K. From 
experience with YNi2B2C we know sharp cusps in (7(H) imply point nodes and not line nodes |l05j . Therefore the 
specific heat for T < 0.3K suggests the appearance of point hke minigaps for k || (1, 0, 0), (0, 1, 0) etc. These may be 
reproduced with the secondary energy gap 

A2(k) ~ e=^*"*cosx(l - acos(4(/)) cos^ x) 
with a < 1. Then the secondary gap has the same symmetry as a 2D /-wave superconductor. 

Exercises 4. 

4.1. Calculate g(H) when H is rotated within the z — x plane. Consider here / = cosx, sinx and cos(2x)- 

4.2. Calculate g(H) when H is rotated within the x — y plane. Consider the same set of /'s as above. 



VI. MAGNETOTHERMAL CONDUCTIVITY 



Experiments on the angular dependent thermal conductivity in YBCO [106111071110811091 1 indicate that the thermal 
conductivity is sensitive to the nodal directions, although theoretical interpretation of this data has been given only 
recently 65]. In order to calculate the thermal conductivity it is crucial to incorporate the effects of impurity scattering 
and of the Doppler shift on an equal footing. Then there arise natually two limiting cases: a) the superclean limit 
(rA)i/2 « |v • q| > and b) the clean limit {TAy/^ X |v • q| » F. 

First it is necessary to determine the quasiparticlc lifetime in the presence of impurity scattering and the Doppler 
shift. This is given by |^] 

ttF 2 

Co = -r[<Coln(-===) + :rtan-i(a;/Co)>]-i (53) 
2A V<^o + ^ 

where ACq = ImiCo) at a; = is the quasi-particle relaxation rate on the Fermi surface and x = |v • q|/A. This 
formula applies to the class of f's introduced in II, which include d-wave superconductivity with / = cos(20). Also 
in the later analysis we assume that H lies in the ab-plane. 



A. Superclean limit 



In the superclean limit we can assume that a; ^ Cq. Then Eq.(52) is solved as 



Co = ^<x>-' --{^f < X [ln(2/x) - 1] (54) 
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On the other hand, in the clean Hmit we obtain 



ttT 1 



Also when H is perpendicular to the a-b plane we obtain 

<.^>^f^ln^ (56) 

The thermal conductivity within the a-b plane and in the limit T <C A is given by 

/ 1 + ^^gg^^';/-^.^^ \ 

K\\/l^n{= Kxx/Hn = Kyy / Kn) = — ( — , , ■ \9 , j? / (5''') 

4A \ Re^/(Co + ixy ^ r I 

where k„ — ""^^m thermal conductivity in the normal state. Then in the superclean limit Eq.(56) reduces to 



< x>= - <x>'^ (59) 



4ACo 4 
1 v^{eH) 
^ A2 



(60) 



In the last step we assumed H || c. The H-linear thermal conductivity was first observed in Sr2Ru04 at T = O.SivT |67j . 
which indicates nodal superconductivity with / — e^"^ cos(x) in Sr2Ru04 and that the system is in the superclean 
limit. More recently the H-linear thermal conductivity was observed in PrOs4Sbi2 'ti'I for T < 0.3K. 
Now in the presence of a magnetic field in the a-b plane, we must first generalize Eq. 57 as |63| 

Kxx/ku — — < X >< (1 cos(20 ))a; > (61) 

TT / 
Kyy/Kn — — < X >< (1 — cos(2(/) ))x > (62) 



and 



TT / 

Kxy/nn = ^ ^ '^'^ sin(20 )x > (63) 



where the angle <f) refers to the direction of the quasiparticle wave vector within the a-b plane. Then for / — cos(20) 
Eq.(58) gives 

v^eH 

n\\/nn{= Kxx/ku = Kyy/K„) = (0.955 0.0285 cos(4(/.))2 (64) 

and 

TT^A 



Hxy/nn = -33-2(0.255 sin(20)(O.955 + -1-0.0285 cos(4(/))) (65) 



Here v = ^Jw^- Theref ore k \\ should exhibit the fourfold term with a magnitude variation of the order of ^ 6%. 

In early experiments '106*107'109] a fourfold term of comparative magnitude was found in YBCO but of opposite 
sign. Similarly Ocaiia and Esq uinazi ,1Q8J found the Hall thermal conductivity ^ sin(20), but of opposite sign. These 
problems were clarified in |63- AH experiments on YBCO were performed at T |v • q| >, whereas the present 
theory appfies only for T <C< |v • q| >. 

Also, thermal conductivity measurements in CeCoIns 68] and k-(ET)2Cu(NCS)2 [i^ have seen evidence for a 
fourfold term. If we assume that these experiments are done in the region T <C< |v • q| >, we have to interpret 
the data as dxy-^SLve superconductivity and dx^-y^ symmetry for CeCoIns and k-(ET)2Cu(NCS)2 respect ively . The 
former identification appears to be consistent with the magnetospecific heat data of CeCoIns reported in [llOj . For 
/ = e^"^ cos X, we obtain [6^ 



Kxx/Kn = ^2 X 1.479(1 -0.0416 cos(20)) (66) 



The magnitude of the twofold term seen in [63 is consistent with Eq.(66). 
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B. Clean limit 



In the clean limit Eq.(57) is transformed as 



K||/ko 



< > 



67r2 



■In 2a/ — 




2A \ v'^{eH) 

rA 



ln( 



(67) 
(68) 



where kq = i^\\{H 0) and we have assumed H || c. In the clean limit the field dependence is given by /v|| 
HIti{Ho/H). In a magnetic field within the a-b plane we obtain 



Kxx/kq = 1 + In I 2a 



/2A\ v^(eH), , A , , 

' ^ -[ln( ^ ) -0.072 + 0.041 cos(40)] 



rA 



vVeH 



(69) 



and 



Kxy/Ko = —2 < sin(2(/))a;^ > 
6Cq 



-^4— ln(2A/ — )sin(2</.)[ln(^=)-0.42] 
Stt^I a V ttI wVei/ 



(70) 
(71) 



where we have assumed / = cos(20). 

The b-axis thermal conductivity of K-(ET)2Cu(NCS)2 determined by Izawa et al jH^ exhibits not only the four-fold 
term predicted in Eq. 65, but also a two-fold term that was not predicted. We can interpret the two-fold term 
as being due to the admixture of an s-wave component. The experimental data is most naturally interpreted as 
A(k) ~ cos(20) - 0.067 JUj. Finally for T >< |v • q| > we obtain 65] 



10 ^A^ ^ n' ^1.76T" 



(2^) 



2A i^eH) 
^I.IQT' A2 



ln(^i^--^(l-cos(40)) 
v\/eH 16 



1 ■ M^'^(e^) , / 2A , 4A , 

Kxy/Hn = ^ Sm(20) IQIt^) M-T^^) 



vy eH 



(72) 
(73) 



It appears that bot h Eq. (67) and (68) describe consistently the angle-dependent magnetothermal conductivity observed 

in YBCO lioalioal . 

Exercises 5. 

5.1 Indicate how to derive Eq.(58) from Eq.(57). 

5.2 Consider k^z for f's with horizontal nodes (i.e. / = sinx, cosx, and cos2x) when the magnetic field is rotated 
within the z-x plane. 



VII. SUPERCONDUCTIVITY WITH POINT NODES 



A. Borocarbides YNi2B2C and LuNi2B2C 

This class of superconductivity was discovered in 1994 |112| . Their relatively high transition temperatures (Tc — 15.5 
K and 16.5 K for YNi2B2C and LuNi2B2C, respectively) and their interplay between magnetism and superconductivity 
are of great interest. In the following we will focus on the two superconducting borocarbides YNi2B2C and LuNi2B2C, 
which have no indication of magnetism. 

The presence of a substantial s-wave component in the superconducting order parameter A(k) in YNi2B2C was 
established by subst ituting Ni by a small amount of Pt. An opening of the energy gap was observed by specific heat 
measurements |ll3j . On the other hand, this superconductivity exhibits a number of peculiarities unexpected for 
s-wave superconductors. For example the \fH dependence of the specific heat in the vortex state |ll4lll5j | indicates 
a nodal superconductor. Furthermore the presence of de-Haas-van-Alphen (dHvA) oscillation in the vortex state in 
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LuNi2B2C down to H = 0.2Hc2 suggests again nod al su perconductivitv. |ll6lll7| Also T-^^ in NMR exhibits a 
power law consistent with nodal superconductivity |ll8l |. In addition A(k) exhibits a fourfold symmetry as seen 
froni the angular dependence of the upper critical field Hc2 when the magnetic field is rotated within the a-b plane 
[n? 120]. Also, unlike the superconductors discussed in previous sections, here we are dealing with a superconductor 
with a 3D Fermi surf ace. 
Then we postulate |105| 

A(k) = (A/2)(l - sm\e) coam)) (74) 

where 6 and (f) are polar coordinates describing k, the quasiparticle wave vector. This order parameter belongs to 
the hybrid representation of s+g-wave superconductivity. The precise matching of the s and g-wave components is 
necessary in order to find the nodal excitations as observed experimentally. Within a model Hamiltonian the stability 
of such a precise matching is considered in [12LJ . The corresponding A(k) is shown in Fig. 13. The quasiparticle 
density of states is given by 



FIG. 13: Order parameter for the borocarbides YNi2B2C and LuNi2B2C 

G{E)^\x\ f dyF{y)Re ^ \ (75) 
Jo \x — V I ' 

where x = E / and 

1 dz 

with uq ^(1 - (|1 - y|)i/2). We note that F{2 - y) = F{y). The DOS is shown in Fig. 14. For \E\/A = a; < 1 
the quasiparticle density of states can be approximated by 




E/A 



FIG. 14: Quasi-particle density of states for the borocarbides YNi2B2C and LuNi2B2C 
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9{E) = ^il + i9/i8nm\E\/A) + ...) 



which gives 



27C(3) 



(T/A) 



where v — ^VaVc and 



i((l - sin^ eisin^ + (1 - sin^ Q 



cos 



(77) 



(78) 



Unhke the point nodes discussed in Ref. T|, the point nodes here are quadratic which gives N{E) ^ \E\ and Cg ^ 
In the presence of a magnetic field the quasiparticle density of states is given by 



(79) 



(80) 



Here 9 and cj) are the direction of the magnetic field where the polar axis is taken parallel to the c axis. We show in 
Fig. 15 1(6, (j)). In particular for 6 = tt/2 the DOS exhibits cusps at (p = 0, 7r/2, tt, 37r/2. These cusps appear also in 
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FIG. 15: The angular function I{e, 



the specific heat, since 



Cs/jnT = g(0,77) 



(81) 



for T ^ A. Indeed, such cusps have been seen by Park et al |l22l ] Unfortunately, however, their data are limited to 
6 = tt/2 and cannot explore the whole extent of 1{6, (j)). 

Due to the presence of an s-wave component, the effect of impurity scattering is very different from the usual nodal 
superconductors discussed in III and IV. (123i | First of all there is no resonant scattering. So the Born approximation 
suffices. Also, unlike in other nodal superconductors the energy gap opens up immediately in the presence of impurities. 
The energy gap ujg is given in a good approximation bywg =:r(l + -£)~^. Therefore the specific heat and the thermal 
conductivity decrease as T ^ like Cs/T ~ k/T ~ {ujg/Tye^^^^'^ . Hence, there is no universal heat conduction. 

The gap equation is solved in the presence of impurities. We find 



ln(T./T.o) = ^[*(- 



r 1 
-*(-)! 



(82) 



which is compared with Eq.(37) for usual nodal superconductors. We show in Fig. 16 Tc/Tco, A(0, r)/Aoo versus 
F/Aqo. Both Tc/Tco and, A(0,r)/Aoo decrease much more slowly as F/Aqo increases. 
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FIG. 16: Tc and A(0,r) for borocarbides with impurities 
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FIG. 17: Thermal conductivity in vortex state in YNi2B2 



Finally the thermal conductivity in the vortex state is given by |12 



- y')'/' ~ ^j;(arccosy - y{l - yY^'Ml - y) 



2\n{2/x) 



(83) 
(84) 



where x — — — ' ^ ,y = We show in Fig. 17 data by Izawa at al |69|, where the magnetic field is rotated 
conically around the c-axis. The clear cusps for 6 — tt/2 disappear as 9 is decreased. This is the clear sign of point 
nodes at k = (100), (010), (-100) and (0 - 10). The effect of Pt-substitution of Ni in Y(Nii_^Pt3;)2B2C with x=0.05 
is studied by Kamata et al '12^. As seen from Fig. 19 the angular dependence in Kzz disappeared completely for 5% 
substitution of Pt. From Eq.(78) and Tc = li.lK of the 5 percent Pt substituted system, we can deduce F = 23.8i4r. 
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FIG. 18: Effect of impurities on thermal conductivity - theoretical 
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FIG. 19: Effect of impurities on thermal conductivity - experimental 



This r is much larger than Tco- Clearly parallel experiments with Pt-substitution less than 1% arc highly desirable. 
The present A(k) describes very well the ultrasonic attenuation data of YNi2B2C ^126, 122,]. 



B. Skutterudite PrOs4Sbi2 

T he skutterudite PrOs4Sbi2 is a heavy-fermion superconductor with transition temperature Tc ^ 1.8K |l28ll29l 
ll3Cll |. The angle-dependent magnetothermal conductivity data of this material inidca tes a n interesting multiphase 
structure characterized by the gap functions AA(k) and AB(k) with point nodes. |70ll7llll3lj More recently, there has 
been mounting evidence for triplet superconductivity in this compound. First, from /x-SR measurement Aoki et al 
|132| discovered a remnant magnetization in the B phase of this compound indicating trip let pairing. Also more 
recent thermal conductivity data by Izawa et al is not c onsistent with the singlet model 131], but is consistent with 
p+h-wave superconductors tTlj]. Finally, Tou et al |l33l | reported NMR data of PrOs4Sbi2 of which the Knight shift 
suggests the triplet pairing. 
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The triplet pairing together with the position of point nodes give almost uniquely |7l| 

AA(k) = ^dAe±*^'(l-k/-k2^~k3") (85) 
AB(k) = dAe±^"*^(l~k3*) (86) 

where e^*"^! = (/fcs ± ik3) / y^{kf+kl) , etc. Here (p, in Eq.(85) is one of (/>! , </)2 and 03 . Also we have chosen the 
nodal direction of Ab parallel to (0 1) in Eq.(86). Also |AA(k)| has the cubic symmetry consistent with the 
crystal symmetry of PrOs4Sbi2, while |AB(k| has the axial symmetry. We show in Fig. 20 |A^_(k| and |AB(k)|. 
Very recently the magnetic penetration depth of the B-phase in PrOs4Sbi2 has been reported |l34j . They applied 





FIG. 20: Proposed A-phase (left) and B-phase order parameters for PrOs4Sbi2 

a magnetic field parallel to each of the 3 crystal axes and determined the magnetic penetration depth and found 
isotropic supcrfluid density. At first sight this is clearly in contradiction to our AB(k) in Eq. (86). However, this is 
understood, if we assume that the stationary magnetic field controls the symmetry axis of AB(k) |l35j . Indeed the 
ground state energy is favorable when the nodes are aligned parallel to H. Also the dependence of the superfluid 
density for T <^ A follows from this assumption. As in all the triplet superconductors, our superconducting order 
parameter breaks the chiral symmetry. So both order parameters are six-fold degenerate. 
The quasiparticle DOS of these states is given by 

for X — This is evaluated numerically and shown in Fig. 21. For \x\ ^ 1 we obtain 

g{E) ~ Tr\x\/A A-phase (88) 
~ 7r|a;|/8 B-phase (89) 

In the B-phase there is a logarithmic singularity at = A as in d-wave superconductors, whereas in the A-phase the 
singularity is split into 2 cusps at E/A = 3/4 and 1. In the vortex state the quasiparticle DOS is given by 

1 v\feH 

5^(H) = 7^(0,0) A-phase (90) 

1 v\feH 

-lB{0,(j)) B-phase (91) 



4 A 



where 



/a(6I,0) = sin^-l- (1 -cos^ 6* sin^c/))^/^ + (1 -cos^ 61 cos2</>)^/2 (92) 
Is (61,0) = sin 6*. (93) 
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FIG. 21: Quasiparticle density of states for p+h-wave superconductor 



In the B-phase we have assumed that the nodes are parallel to (0 1) 
are given by 



Then the low-temperature specific heat, etc. 



CsHnT = 5(H), xs/XAf = .g(H) 
Ps(H)/ps(0) = l-.g(H) (A-phase) 
Ps||(H)/ps||(0) - l-3ff(H) (B-phase) 



(94) 
(95) 
(96) 



The last expression means the superfluid density measured parallel to the nodal directions. 

In p+h-wave superconductors the impurity scattering is handled similarly to other nodal superconductors. For 
example, the superconductivity in the A-phase exhibits the universal heat conduction with kqq/T = 7r^n£'i?/12A(0). 
On the other hand in the B-phase the heat current has to be parallel to the nodal direction. Then we will have 
i^oo/T — TT^nEp/8A{0). When the heat current is perpendicular to the nodal direction in the B-phase, the thermal 
conductivity vanishes like T^. 

Now let us consider n^z in the vortex state when H is rotated in the z-x plane with 9 the angle H makes from the 
z axis. Then in the superclean limit ((FA)^/^ <C v^/eH) we obtain 



sm A — phase 



(97) 



64A2 



sin^ 6 B — phase 



where in the B-phase we assumed that the nodes are parallel to (0 1). Similarly in the clean limit ((FA)^/^ ^ v\ 
we obtain 



Kzz/k.Q 




12FA 



I2A 

— )sin^6lln(- 
F v\ 



eH sin 9 



) (A — phase) 
(B — phase) 



^) 

(99) 
(100) 



where kq = Kzz{H = 0). In Fig. 22 we show the experimental data at T = 0.35 K together with the theoretical 
fit from Eq.(90) and (91). As is readily seen the fits are excellent with A/v^/eH — 5 and 3 for data at H = 0.5 T 
(B phase) and 11= 1.2 T (A-phase) respectively. Then making use of the weak-coupling theory gaps A^(0) — A.2K 
and Ab{0) = 3.5K for the A-phase and B-phase respectively, we can deduce v = 0.96 x 10^ cm/sec and F ~ O.IK. 
These values are very reasonable. From the de Haas- van Alphen measurement v is estimated to be 0.7 x 10^ cm/sec 
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FIG. 22: Thermal conductivity of PrOs4Sbi2 in the vortex state 



(a band), 0.66 x 10^ cm/sec (/? band), and 0.23 xlO^ cm/sec (7 band) |136| . Perhaps the triplet superconductivity 
in PrOs4Sbi2 is n ot so surprising. The superconductivity in Sr2Ru04, UPts, UBeis, URu2Si2 and UNi2Al3 appears 
to be triplet |l37j . More surprising is the high degeneracy of the ground state due to the chiral symmetry-breaking. 
Both the A-phase and the B-phase have six-fold degeneracy. This means that a variety of topological defects in both 
the A-phase and the B-phase are likely, which deserve further study. The multi-phase structure is rare but at least 
we have seen it in superfluid '^He and in UPta. 

On the o ther hand the exact location of the A-B phase boundary is still controversial. A recent therniodynamical 
study |l38j| indicates that the A-B phase boundary is rather parallel to the upper critical field Hc2 of the A phase. 
Somewhat surprisingly the less symmetric B phase is realized for the low-field region. Is this related to the antiferro- 
magnetic quadrupolar (AFQ) state which appears above iJ ~ 4T ^139j ? It looks like both the superconductivity in 
YNi2B2C and in PrOs4Sbi2 have great futures. 



C. Summary on A(k) 



a). Quasi-2D systems 

1. d-wave superconductors (/ = cos(2(?!)); hole doped high-T^ cuprates (YBCO, LSCO, Bi2212, T12201), electron- 
doped cuprates (NCCO, PCCO), CeCoIns and k — {ET)2Cu{NCS)2- As to the gap symmetry of CeCoIns there is 
controversy between d3,2_j,2-wave versus d^^y-wave. 

2. Sr2Ru04. Early experiments established triplet pairing and chiral symmetry breaking |l4Cil4llll42j |. Whe n the 
high quality single crystals became available, both the specific heat data 58] and the magnetic penetra tion d ata |143| | 
exhibited characteristics of nodal superconductors. Therefore the p-wave superconductivity propo sed i n |145| is clearly 
out. In order to save this situation amultigap model was proposed by Zhitomirsky and Rice |l4fj |. However, our 
analysis of the optical conductivity |l47j indicates that there is little room for p-wave superconductivity. Therefore 
the important question is: where are the line nodes in A(k). 

Clearly the magnetothermal conductivity data (66] and the ultrasonic attenuation data |l48j ] support horizontal 
nodes as in / = e^*"^ cos x- Therefore, except for the specific heat data by Deguchi et al |l04j we mentioned earl ier al l 
available data are consistent with the 2D f-wave model. Also we have proposed that the optical conductivity jl47j . 
the Raman spectra |l49l ]. and the supercurrent experiment jl5^il51jlbelow T =0.1 K will provide the definitive test 
of A(k) in Sr2Ru04. The review paper by Mackenzie and Maeno p52l ] is excellent, but does not cite the significant 
evidence for nodal superconductivity in Sr2Ru04, if we concentrate on the data from high quality single crystals. On 
the other hand, we agree that definitive experiments below T = 100 mK are highly desirable. 
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3. UPd2Al3 As briefly described in IV, we have proposed f = cosf2Yl.|89l| 

4. We are awaiting the magnetothermal conductivity data of UNi2Al3, URu2Si2 and CeCu2Si2. 
5) 3D systems 

a) . YNi2B2C, LuNi2B2C. A(k) = (A/2)(l - sin^ 6 cos{A(j))) 

b) . PrOs4Sbi2, p+h-wave superconductivity. 

VIII. D-WAVE DENSITY WAVE 

In the preceding lectures we have reviewed a variety of nodal superconductors which have appeared since 1979. 
For some of them we have succeeded in identifying the gap symmetry. All these analyses are based on the BCS 
theory of nodal superconductors. Also the semiclassical approximation played an important part in elucidating the 
quasiparticle properties in the vortex state. 

Th erefore it is very natural to contemplate a parallel development in charge density wave (CDW) and spin density 
wave |2ll53) systems The unconventional density wave (UDW) was speculated as a possible electronic ground state 
in excitonic insulators in 1968 154J. However, the recent surge of interest in UDW is in part due to the pro posal 
that the pseudogap phase in high-Tc cuprate superconductors is a d-wave density wave fd-DW) . |4^ll4ll45ll4a | It is 
important to point out that two kinds of d-DW have been considered. The first one is commensurate and individual 
square lattice enclosing a circulating current 43 45j analogous to the flux phase Ji]. Therefore the ground state has 
Z2 symmetry and "visons" as topological defects. On the other hand we consider a d-DW, which is in general 
incommensurate. [46) Therefore as in a superconductor the ground state has U(l) symmetry associated with the axial 
gauge transformation 

- e-^/^cI^Q (101) 

The d-DW can have "phase vortices" similar to the ones in conventional CDW [i^ . 

Here are a few characteristics of unconventional density waves (UDW) |l55j . First of all, the transition from the 
normal state to UDW is a metal-to-metal transition. Though the quasiparticle density decreases due to the opening 
of a partial energy gap, UDW are conductors down to T = K. In high-Tc cuprates the resistivity ~ in the 
normal state changes to the resistivity ~ T for example. Also, we shall see that the quasiparticles in UDW are a 
standard Fermi liquid. Therefore this anomalous resistivity behavior does not imply the existence of a "non-Fermi 
liquid". Since (A(k))=0, there will be no x-ray or neutron signal for charge density or spin density. Therefore UDW 
is sometimes called a condensate with a hidden order parameter |45). 

The quasiparticle Green function in UDW is given by 

g-\k,uj) ^ Lu - r]{k) - ^' {k)p:i ~ Aik)pi (102) 

where ^ (k) — |(^(fc) — £,{k — Q)),rj{k) — \{S,{k) + £,{k + qj) and Q is the nesting vector. Also the p's are the Paufi 

matrices operating on the spinor space made of |c^) and \c\_q). Here we consider only UCDW for simplicity and 

spin indices are dropped. In the following we drop the prime for ^ (k), since there will be no confusion. Then the 
quasi-particle energy is given by 

uj = T](k) ± v/|^(fc)TA2(fc) (103) 

which is identical to that for nodal superconductors except for r]{k), the imperfect nesting term. 
Here are two urgent questions: 

Where can we find UDW? If they exist, how can wc identify them? As to the second ques tion, we believe that two 
hallmarks of UDW are the angle-dependent magnetoresistance and the giant Nernst effect [iSiij . Bo th of the se are 
consequences of the Landau quantization of the quasiparticle spectrum in UDW or Nersesyan's effect. |T5^157| 

A. The Nersesyan effect 

In 1989 Nersesyan et al pointed out that the quasiparticle motion in UDW is quantized in the presence of a magnetic 
field perpendicular to the conducting plane. Then the QP spectrum becomes 



^'q± ^y2nvv2e\Bcos0\ (104) 
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where n=0,l,2,3... and the n ^ states are doubly degenerate and V2/V = A/Ep |l58j |. Here A is the maximal energy 
gap in UDW and we assumed 77 (the chemical potential) independent of k. First the magnetoresistance is given by 

RiB, or = 2a. + 4., ( ^-^^^^^) (105) 

\coshxi + cosh^o / 

where Xi = (3y^2evv2\B cos 6*1, ^0 = and /3 = l/ksT. Here we considered only the 2 lowest Landau level s. A 
slight modification of this simple formula can desc ribe the low-temperature phase of a-(ET)2KHg(SCN)4 jl59l | and 
the metallic phase in the Bechgaard salts |l6Cil6ll| . See Fig. 23. Second, in the presence of an electric field the QP 




T(K) p 

FIG. 23: Phase diagrams for low-temperature phase of a-(ET)2KHg(SCN)4 |l5fl|| ('left'l and the metallic phase in the Bechgaard 
salts (right) 



orbits drift with vd — (ExB)/|i3p. This gives rise to a large negative Nernst effect jl58l |. We obtain 

SR 

a^y = = 2ei?(ln2 + 21n(2cosh(xi/2))-.Titanh(a:i/2)) (106) 

B 

where S is the en tropy associated with the QPs. Indeed a large negative Ner nst effect has been observed i n a- 
(ET)2KHg(SCN)4 Te^, in the pscudogap region of high-T^ cuprates |l63ll64llT65| . in the CDW region of NbSca Tel 
and more recently in CeCoIus 167 168). Recently we have analyzed the Nernst coefficient observed in CeCoIns 163| 
in terms of Eq. 105. This is shown in Fig 24 a) and b). We obtain an excellent fit as is seen readily. On the other 
hand, in order to describe the temperature dependence we have to assume 

A{T)v{T) = a + bT^ (107) 

whose origin is unclear. However, it is possible that the above temperature dependence indicates the presence of a 
quantum critical point (QCP) in CeCoIns. 

In add ition, we have shown recently the giant Ne rnst effect observed in the underdoped region of YBCO, LSCO and 
Bi-2212 '163, 165] are described in terms of dDW 16!^. Very recently a similar but positive giant Nernst effect has 
been reported in the AF phase in URu2Si2 .IJQl. Inde ed there has been a suggestion that CDW in NbSe2 is UCDW 
[ml and the AF phase in URu2Si2 is USDW |l72iT73l |. So far we have identified 7 candidates for UDW which include 
the 3 high-Tc cuprate superconductors YBCO, LSCO and Bi-2212 in the underdoped region. 

Coming back to high-Tc cuprates, 1) A(k) — Acos(2(/)) in the pscudogap region is determined by ARPES jl74l |. 
2) The similar phase diagram as in Fig. 2 is obt ained when T is replaced by the lo w- tem perature energy gap which 
is measured by STM |l75lll76l |. ARPES jl4ll77j and the universal heat conduction flTSj. In particular, if one puts 
A = 2.14T* these two phase diagrams match almost p erfectly . Needless to say, A/Tc = 2.14 is also valid for d-wave 
density wave in the weak-coupling limit when |ry| ^ A |l7l74j . 



IX. GOSSAMER SUPERCONDUCTIVITY 



In the first lecture we discussed Laughlin's critique 51] on Anderson's RVB state (Eq.(2)). The new wave function 
suggested by Laughlin can be interpreted as d-wave superconductor in the presence of another order parameter '52^531 . 
From the phase diagram for high-Tc cuprates, the most relevent state is a d-wave superconductor in the presence 
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B(T) B(T) 

FIG. 24: The magnetic field dependence of tfie Nernst coefficient is plotted for (left) T = 1.3 K, 1.65 K, 2.5 K, 3.5 K and 4.8 
K from top to bottom, and for (right) T — 7.3 K, 10.5 K and 15 K from bottom to top. The circles denote the experimental 
data, while the solid line is our fit. 



of d-wave density wave. Also this is the prevailing feature in heavy-fermi qn systems li ke URu2Si2, CeCu2Si2 [ij 
CeCoIns, CeRhlns, UPd2Al3, and organic conductors like k — {ET)^ salts |l8Cil8llll82 

Coming back to high-Tc cuprate superconductors, A(k) for d-DW can be readily identified as the order parameter 
phenomenologically introduced by Tallon and Loram |5C| . For example, the superfluid density in the gossamer 
superconductor at T=0 K is given by 



while the c-axis superfluid density is 



A2(0) 

- Ai(0) + AKO) ^'''^ 



/O.c(0,77) - , ^'^"^ = (109) 



where Ai(0) and A2(0) are the maximal gaps of d-wave superconductor and d-wave density wave, respectively, at T= 
K. 

It is of great interest to explore a variety of transport properties in the gossamer superconductivity. For example, 
if T], the chemical potential is negligible {\rj\ <C T), we obtain again the universal heat conduction (Eq.(3)) where now 



A = yAf+A|. Also, ARPES would see the energy gap (17: 

A(k) = ^J (A2 cos(2</>) - r/)2 + A^ cos2(20) (110) 
~ |Acos(20) -?7| (111) 

If you look carefully at the phase diagram of Bi-2212 (see for example |l77j V unlike LSCO and YBCO the supercon- 
ducting dome in Bi-2212 is completely covered by the pseudogap phase. This suggests that the superconductivity in 
Bi-2212 is gossamer for the whole doping range. This may be one of the reasons why A(0)/rc ~ 5 in Bi-2212 is so 
large. In this sense the systematic study of the optical conductivity, the Raman spectra and the thermal conductivity 
(with special attention to the doping dependence of these quantities) is of great interest. 

In summary, in exploring unconventional or nodal superconductors, we have encountered a vast forest inhabited by 
many unconventional density wave (UDW) and gossamer superconductors where these two order parameters coexist. 
Surprisingly, all of these ground states have been expected from the infrared instability of the 2D and 3D Fermi liquid. 
Therefore we can restore the legacy of Landau's Fermi liquid theory and BCS theory of superconductivity in a proper 
persp ective. Then armed with the Green function methodology as in Abrikosov, Gor'kov and Dzyaloshinski (AGD) 
|183l| we will be fully prepared to explore the plethora of new ground states. 
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